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Young's moduli of regular two-dimensional truss-like and eye-shape-like structures are simulated 
by using the finite element method. The structures are the idealizations of soft polymeric materials 
used in the electret applications. In the simulations size of the representative smallest units are 
varied, which changes the dimensions of the cell-walls in the structures. A power-law expression 
with a quadratic as the exponential term is proposed for the effective Young's moduli of the systems 
as a function of the solid volume fraction. The data is divided into three regions with respect to the 
volume fraction; low, intermediate and high concentrations. The parameters of the proposed power- 
law expression in each region are later represented as a function of the structural parameters, unit- 
cell dimensions. The presented expression can be used to predict structure/property relationship in 
materials with similar cellular structures. It is observed that the structures with volume fractions 
of solid higher than 0.15 exhibit the importance of the cell- wall thickness contribution in the elastic 
properties. The cell-wall thickness is the most significant factor to predict the effective Young's 
modulus of regular cellular structures at high volume fractions of solid. At lower concentrations of 
solid, eye-like structure yields lower Young's modulus than the truss-like structure with the similar 
anisotropy. Comparison of the numerical results with those of experimental data of poly (propylene) 
show good aggreement regarding the influence of cell-wall thickness on elastic properties of thin 
cellular films. 



I. INTRODUCTION 

In 1930's Munters and Tandbergi invented the first 
man-made cellular materials, which opened a new era 
on high-tech materials. These new class of materials, 
called foams, were introduced in heat insulation and 
packaging products in early days. Nowadays, they are 
widely used in variety of engineering applications 2 ^, such 
as aeronautics, impact absorbing materials^, electro- 
mechanical sensors&l, etc. Recently, cellular polymers 
with electrical insulating properties have been shown to 
exhibit electro-mechanical activitjiSiSiSiiSiiiiiS and are 
called "ferro-electrets" . Amplitude of stored charge in- 
side the foam and foam's elastic property are the main 
parameters that influence their electro-mechanical prop- 
erties. The structure of these cellular materials is lay- 
ered polymer sheet with voids (resembles a puff pas- 
try) , scanning electron micro-graphs have been presented 
elsewhereZi&^ii&ii. The materials are produced as films 
with lateral dimensions much larger than the thickness 
one. There are possibilities to alter their internal struc- 
ture with post pressure processing^ by changing the di- 
mensions of voids. Therefore, in this paper we focus on 
the mechanical properties of cellular structures use in the 
ferro-electret applications, and assume two regular ge- 
ometries, which resemble the structure of these cellular 
polymers. 

The elastic properties in the thickness direction, ef- 
fective Young's moduli, of two these structures are cal- 
culated using the finite element (fe) methodiS*iSiiLi2iiS 
and compared with each other. The obtained effective 
Young's moduli of the structures are expressed by a 
power-law function similar to that proposed by Gibson 
and Ashby—, however, here the exponent is expressed as 
a quadratic function of solid volume fraction q, (q = p* / p 



with p* and p s being the densities of the cellular and solid 
materials). The elastic properties of the structures are 
expressed in the whole volume fraction range with the 
proposed expression as a function of structural param- 
eters. Experimental data from literature are also pre- 
sented for comparison with the numerical model. 

The cellular structure is a solid-gas mixture, a hetero- 
geneous medium, whose macroscopic or effective proper- 
ties can be expressed as a function of properties of the 
constituents and their volume fraction, 

<P e = (1) 

where is the effective property sought and subscript 
'i' represents the properties the constituents with volume 
fractions qi, = 1. For voided solids, to predict the 
effective Young's modulus E e is nontrivial, due to dif- 
ficulties in assigning the mechanical properties for the 
void. However, in such cases numerical simulations are 
standard methods. When the voids are open, meaning 
that they do not confine incompressible fluids, we can ne- 
glect this phase in the calculations. Gibson and Ashby 2 
have used beam theory to obtain the famous power-law 
relation between the effective Young's modulus E e and 
the volume fraction of the solid material q 

E e = <Lq n E s (2) 

where, £ and n are parameters depend on the micro- 
structure, and E s is the Young's modulus of the solid. 
The value of n lies between 1 and 4, yielding a wide range 
of effective properties for a given solid volume fraction 
<j 20 . Experimental observations on various open-cell cel- 
lular materials suggest n = A Since the dependence of 
effective properties on micro-structure are not well under- 
stood in composite materials, the exact form of £ and n 
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FIG. 1: (a) Truss-like and (b) eye-like structures. The shaded 
regions are the representative units used in the calculations 
with the appropriate boundary conditions as shown in Fig. |2] 



are not known. Therefore to optimize and to predict foam 
properties numerical simulations and structure/property 
investigations play an important role in the current re- 
search activities2*i2i2i. 



II. CONSIDERED GEOMETRIES 

In order to simulate the elastic properties of the 
cellular structures in ferro-electrets presented in the 
Iiterature£*2i2iiii4 two geometrical models composed of 
repeating unit cells on a regular triangular lattice are 
constructed. In the current study influence of the struc- 
tural differences are investigated. In the model geome- 
tries, the solid- void systems consist of straight and curved 
layer structures as presented in Fig. ^ The layers are 
connected at the triangular lattice points. In this model, 
the voids become lozenge and eye-shaped geometries, as 
illustrated in Fig.QJi and H 3 , for truss- and eye- like struc- 
tures, respectively. The structures are refered as truss- 
like and eye-like in the text. The marked rectangular 
regions [(0, 0) < (x, y) < (a/2, 6/2)] in Fig.UJi andHt are 
the unit cells (or the computational domain f2) used in 
the FE calculations, 

V = (0, 0) < (x, y) < (a/2, 6/2) A (x, y) € U s . (3) 

The parameters a and 6 are used to model different void 
geometries-cellular structures. 

The cell-wall thickness of the truss-like structure t* is 
calculated from the dimensions of the unit-cell a and 6 
and volume fraction of the solid q. 

t*(g,a,6) = 2V34sin(2a) (4) 

with 

A = (6/2 - to tana) 2 + (a/2 - to) 2 
to = \J\ — qab/(4 tan a) 
a = arctan(6/a) (5) 



FIG. 2: Unit cell (representative volume element) in the sim- 
ulations. The region inside the dashed line (- -) is con- 
sidered to be mirrored in all directions. The dark-gray re- 
gion is the solid medium fl 3 . The void is presented with 
f2„. The boundary conditions are also presented as a bound- 
ary force dfl* and constrained displacements dfl u with arrows 
and triangles, respectively. The free boundaries are presented 
as dn f . The solid region is allowed to deform only in the 
y-direction. The dimension of the unit-cell is a x b, where 
a = |AB| = |CD| and b = |AD| = |BC|. 



The cell- wall thickness of the eye-like structure is mod- 
eled with the following line expression with wall thickness 
t e : 

k = 6/4[l + cos(27re/a)] + (-1)^*72 

with i = {l,2}, (6) 

where l\ and I2 are the lower and upper lines, respectively, 
which represent the boundaries between the void and the 
cell- wall. The region occupied by the solid, f2 s (x, y) G O, 
is calculated with the help of the intersection points of 
line h with y ~ 0, and of line I2 with y = 6/2, labeled as 

Xa , 

x Ci = a/2[(i-l)- l/7rarccos(2i7&- 1)] 

where i = {l,2} (7) 

The solid region is described by the discrete points (x, y) 
that satisfy the following condition: 

n s = [(x, y) < i 2 ] n l(x, y) > h] n [{x, y) e fi] (8) 

In order to calculate the concentration (or relative den- 
sity) q, the area of the solid region, Q s , is determined 
numerically with an fe procedure, in which the region 
£l s is divided into triangular subregions. The area of f2 s 
is iterated with selected thicknesses t e , until the desired 
concentration q is obtained. 



III. THE NUMERICAL ANALYSIS 

Our fe method is based on the minimization of the 
elastic energy in conjunction with irregular meshing 15,22 . 
The Newton's second law in the two-dimensions, which 
is in reality a system of two linear equations known as 
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TABLE I: Boundary conditions in the fe. The constrained 
edges have u = or v — conditions for 'yes' labels. All 
the other edges are 'free' to move in corresponding directions. 
The applied load per unit length is presented by Fq. 

dtf dn u dfl v dsW 
Load F 

x-constraint yes yes free free 

y-constraint yes free yes free 



us to have small discretization units even for very thin 
cell-walls thicknesses efficiently. In the simulations, the 
volume fraction of the solid and the dimensions of the 
unit-cell are varied. 



IV. RESULTS AND DISCUSSION 



A. Effective properties 



Navier's equations, are solved for spatial points obtained 
from the fe triangulation of Q s . 



V T = K 



(9) 



Here T is the mechanical stress tensor and K is the load 
vector. The stress tensor has three components in two- 
dimensions, which are normal stress components T x and 



and the shear stress component S 



dx 

dTy 

dx 



dS 



xy 



dy 

as 



xy 



dy 



(10) 
(11) 



Including the Hooke's law between the stress and the 
elastic strain, T = cVu, lead to a partial differential 
equation expressed in the global displacements u(u,v). 
Then, Eq. JUJ becomes as follows 



V • cVu = K, 



(12) 



where the coefficient c is a function of the solid mate- 
rials Young's modulus and Poisson's ratio. Due to the 
dimensions of ferro-electret films, the contribution of the 
Poison-ratio is modest to elastic properties in the lateral 
dimensions. Our experience shows that the displacement 
in the thickness mode is much larger than the laterial 
ones. Therefore, the boundary conditions are selected 
such that the structures are not allowed to move in the 
lateral dimensions. In Fig. [3 general geometrical consid- 
erations together with boundary conditions for the nu- 
merical model are presented. The Young's modulus of 
the structure is calculated from the average deformation 
of the edge dfl 1 in the y-direction J gQt v dx. 



E e = bF / vdx 



(13) 



where, b is the total thickness of the unit-cell in the y- 
direction. The solid region fl s in Fig. [5] is used in the 
simulations together with the boundary conditions pre- 
sented in Table |U Meshing of the computational domain 
is an important factor in the fe analysis, we have there- 
fore adapted a procedure which used the cell-wall thick- 
ness t as the meshing parameter. In the calculations, the 
fe triangle sides are adjusted not larger than the one- 
tenth of the cell-wall thickness. This procedure allowed 



Simulation results for three different unit-cell dimen- 
sions are presented in Fig. [3] as log-log J3K) and semi- log 
plots for a/b = {2, 4, 8}. Two significant charac- 
teristics are observed (i) as a/b increases, the normalized 
effective Young's moduli of the structures decrease and 
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FIG. 3: Normalized Young's moduli of truss-like (•) and eye- 
like (■) structures; (a) in log-log and (b) in semi-log repre- 
sentations. The dashed lines (- -) indicate the positions of 
volume fractions at which the elastic properties of the struc- 
tures crossover each other. The solid lines ( ) are drawn 

to guide the eyes. 
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(b) at some volume fraction the dependence of the effec- 
tive Young's moduli on q becomes linear in the log-log 
representation. In addition, when q w 0.15 and q « 0.85 
the effective Young's moduli of the truss-like and eye-like 
structures crossover, which are illustrated marked with 
the dashed lines (- -) in the figure. 

The calculated data for ten different a/b ratios are an- 
alyzed with a power-law expression as quadratic form in 
volume fraction q 



E*/E s = €q 



(14) 



where £, a, b and c are parameters dependent on the 
geometrical quantities a/b and q. The data is divided 
into three volume fraction regions; low, intermediate and 
high, which correspond to the regions separated with the 
dashed lines (- -) in Fig.|3Ji and|3j3. The curve-fitting 
analyses indicated that the parameters £, a, b and c can 
be represented as a function of the geometrical parameter 
a/b. We have adopted a polynomial dependence, 



y = a\x 



a 3 x 



a2X 



a\x + ao 



(15) 



where y is anyone of the parameters y = {£, o, b, c} and 
x = — log(a/6). We have obtained very good agreement 
between the results of the polynomial analysis and the 
FE results. The coefficients of the polynomials for each 
parameter in Eq. I|14|) are presented in Table ITU 
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FIG. 4: Dependence of the fitting parameters £, o, b and c on 
a/b ratio. The graphs at the left and right are respectively for 
truss-like and eye-shape like geometries. The symbols (0 — ♦) 
and (o — •) represents truss-like and eye-shaped structures 
respectively. The open symbols are the results of Eq. (1141 . 

and the filled symbols together with the solid lines ( ) are 

the polynomial curve-fitting results. The parameters a, b and 
c are nearly constant in the considered range of a/b ratios 
for the truss-like structure, therefore, the average values are 
presented in these graphs. The coefnents of the polynomials 
are presented in Table ITTI 



B. Low volume fractions of solid 



Intermediate volume fractions of solid 



At low volume fractions of the solid (q < 0.15), the 
truss-like structure yield higher Young's moduli than the 
eye-like one. The fitting parameters are presented with 
open symbols as a function of a/b in Fig.0] In the figure 

the solid lines ( ) represent the curve fitting results 

with the values in Table HTI using Eq. I|15fl . The two struc- 
tures show interesting behavior when dependence of the 
Young's moduli are considered. The truss-like structure 
can be e xpressed by a simp le relation as the one pro- 
posed bv iGibson and Ashbvl Eq. @ with n = c « 1 — 
(aq 2 + b q <C for q < 0.15). In addition, there is a 
linear relation between pre-exponential log €. in Eq. (|14|> 
and log(a/6) for the truss-like structure, see Table ITTI 

The eye-like structure on the other hand do not ex- 
hibit similar behavior as the truss-like structure. The 
pre-exponential can be expressed similar to the truss-like 
structure's one but the constant term is higher. The 
other parameters in the proposed effective formula yield 
quadratic relation to structure parameter log(a/6), see 
Table [n] and Fig. 0] The parameters at the exponential 
term, o, b and c, have minimum or maximum approx- 
imately around a/b = 4 as shown in Fig. 0] Last but 
not least as the concentration of the solid phase is much 
smaller than 0.1, n ~ c, and equal to 1 and ~ 3 for truss- 
and eye-like structures, respectively. 



At intermediate volume fraction of the solid, 0.15 < 
q < 0.85, the eye-like structure has higher elastic mod- 
ulus E* than the truss-like one. Moreover, the eye-like 
structure shows a moderate change with respect to q com- 
pared to the truss-like one, see Fig. |3t>- The truss-like 
structure indicated a bending point (a knee), while the 
eye-like structure does not. The fitting parameters are 
dependent on the unit-cell dimensions a and fo, and cubic 
polynomials can express the dependence as presented in 
Fig. and Table |ITJ 



D. High volume fractions of solid 

As the concentration approach to higher values q > 
0.85, the structures loose their porous character and the 
structural differences are not pronounceable as at the 
other concentration regions. In this region, the truss- 
like structure has higher effective Young's modulus E* 
than the eye-like one. The fitting parameters of the pro- 
posed effective medium equation [Eq. (|14f> ] are shown 
in Fig. H3 The results of the dependency on a/b ratio 
are also plotted in the figure with the polynomial coef- 
ficients presented in Table [n] One very significant dif- 
ference between two structure is the form of the poly- 
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TABLE II: Polynomial dependence of the fitting parameters log £, o, b and c of the power-law expression in Eq. 1141 for three 
volume fraction regions. Below x = — log(a/fo), and y — d4X 4 + a^x 3 + a,2X 2 + a\x + ao, y being one of the fitting parameters, 
y = {£, o, b, c}. 
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FIG. 5: Dependence of the fitting parameters £, a, b and c on 
a/6 ratio. The graphs at the left and right are respectively for 
truss-like and eye-shape like geometries. The symbols (0 — ♦) 
and (o — •) represents truss-like and eye-shaped structures 
respectively. The open symbols are the results of Eq. (f 14p . 

and the filled symbols together with the solid lines ( ) 

are the polynomial curve-fitting results. The coeffients of the 
polynomials are presented in Table IH1 



nomial assumed, cubic polynomials are enough for the 
eye-like structure, while it is forth order polynomials for 
the truss-like one in two of the coefficients. 



E. Cell-wall thickness 

Finally, we have also investigated the influence of cell- 
wall thickness t on the Young's modulus. In order to 
perform the analysis, we isolate the volume fraction q in 
Eq. (0J for the thickness of truss-like structure t* , and use 
the FE modeling to calculate the thickness t e for a given 
volume fraction q for the eye-like structure. We can then 
convert Fig. [3] to constant cell-wall thicknesses for both 
structures. The results are shown in Fig.[7|for log-log and 
semi-log representations. The values are normalized to 
the base length a of the unit-cells by keeping the height 
constant, 6 = 1. It is clear that for some a/b values 
(4 < a/b < 10) and at high and intermediate concentra- 
tions of solid (q > 0.15), both structures elastic moduli 
yield similar values, when constant cell-wall thicknesses 
are considered. The chain lines (- • -) are drawn to in- 
dicate the quadratic dependence, E*(t) = b2q 2 + biq + bo. 
The anisotropy and structural differences do not play 
any significant role at high solid concentrations, and 
the cell-wall thickness is the main factor that influences 
the elastic properties. At lower concentrations of solid 
(q < 0.15), on the other hand, the structural differences 
affect the Young's modulus E* . The influence of cell- wall 
thickness is most significant at high a/b values at low 
concentrations q. There is a minimum for the Young's 
modulus for a selected cell-wall thickness for both struc- 
tures. The position of this minimum is around q = 0.6, 
and it shows a very slight trend to move towards lower 
concentrations with decreasing t/a. It is remarkable that 
the constant cell- wall behavior is very distinct and we do 
not require to divide the g-axis in regions as before. 

In order to illustrate the success of the model, experi- 
mental values from the literature^ are also plotted with 
open symbols in Fig. The data were obtained by al- 
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FIG. 6: Dependence of the fitting parameters £, o, b and c on 
a/6 ratio. The graphs at the left and right are respectively for 
truss-like and eye-shape like geometries. The symbols (0 — ♦) 
and (o — •) represents truss-like and eye-shaped structures 
respectively. The open symbols are the results of Eq. (f 14| l . 

and the filled symbols together with the solid lines ( ) 

are the polynomial curve-fitting results. The parameter £ is 
nearly constant in the considered range of a/b ratios for the 
truss-like structure, therefore, the average value is presented 
in the graph. The coeflients of the polynomials are presented 
in Table |n| 



tering void dimensions of polypropylene samples by post 
pressure treatment. Later the elastic proporties were 
measured using impedance spectroscopy^. The density 
and the Young's modulus of the solid were p s = 1 g/cm 3 
and E s = 1 GPa, respectively. The experimental data 
display a similar behavior as predicted by the fe model- 
ing. Most of the data points sit on the equi-thickness line 
t/a ~ [0.04 — 0.05], which is in good agreement with the 
scanning electron micro-graphs. The proposed power-law 
expression in Eq. I|14fl with the polynomial dependence 
on the structural parameters can be used to prognosti- 
cate the elastic properties of cellular thin films by using 
their scanning electron microscope images. 



V. CONCLUSIONS 

In this paper, we present numerical calculations of the 
elastic properties of two cellular structures which resem- 
ble cellular polymers used in the electret applications. 
The computational domain is constructed from unit-cells 
(repeating units) with solid phase in the form of truss- 
and eye-like shapes. The size of the unit-cells are varied 
in the simulations to introduce anisotropy. The elastic 
properties of the structures are presented for various unit- 




FIG. 7: Normalized Young's moduli of truss-like (•) and eye- 
like (■) structures for various cell-wall thicknesses (a) in log- 
log and (b) semi-log representation. The cell-wall thickness 
are normalized to the base of the unit-cell a. Sixteen thickness 

values are presented, — 2 < log (t/a) < —I. The solid ( ) 

and dashed (- -) lines are drawn to guide the eye. The 
chain line (- • -) is a quadratic relation in q. The open sym- 
bols are the experimental data taken from the literature-*-. 



cell dimensions. The structures' effective Young's moduli 
crossover each other at two concentration levels, q = 0.15 
and q = 0.85. These concentration levels are used in the 
analysis. The most significant difference between the two 
structures is that the considered eye-like structure shows 
much lower Young's moduli than the truss-like one at 
low concentration of solids, which is due to the curved 
cell-wall. Although, at intermediate and high concentra- 
tions, there are visible differences in the Young's modulus 
of the structures when expressed as a function of concen- 
tration, if the cell-wall thicknesses are taken into consid- 
eration, the Young moduli exhibit similar behavior. We 
can therefore conclude that the cell-wall thickness is the 
most significant factor at high concentrations of solid in 
cellular materials. The influence of structural differences 
are only visible at low concentrations and less isotropic 
structures. The calculated effective Young's moduli E* 
can be expressed with a power-law, in which the power 
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exponent is a quadratic expression on the volume frac- 
tion of the solid q. The coefficients of the quadratic term 
and the pre-exponential are expressed as a function of 
the unit-cell dimensions. A full solution for the whole 
volume fraction spectrum is given for both structure, 
which can be of interest for materials researchers deal- 
ing with similar structures as presented ones. Finally, 
computer simulations are excellent tools to better under- 



stand the structure/property relations in cellular layered 
structures. 
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